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ABSTRACT 

Exploiting the presence of kilohertz quasi-periodic oscillations (QPOs) in the timing 
power spectrum, we find that the soft x-ray emission of the neutron-star X-ray binary 
and atoll source 4U 1636-536 modulated at the QPO frequency lags behind that of the 
hard x-ray emission. Emission in the 3.8-6.4 keV band is delayed by 25.0 ±3.3 /is relative 
to the 9.3-69 keV band. The delay increases in magnitude with increasing energy. Our 
results are consistent with those of Vaughan et al. (1997), when the sign is corrected 
(Vaughan et al. 1998), for the atoll source 4U 1608-52. The soft lag could be produced 
by Comptonization of hard photons injected into a cooler electron cloud or by intrinsic 
spectral softening of the emission during each oscillation cycle. 

Subject headings: accretion, accretion disks — stars: individual (4U 1636-536, 4U 1608- 
52) — stars: neutron — X-rays: stars 
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1. Introduction 

Millisecond quasi-periodic oscillations (QPOs) in 
the x-ray emission of x-ray binaries are likely related 
to physical processes occurring near the compact ob- 
ject and may, thus, provide information about accre- 
tion flow in the strong gravitational fields near the 
compact object (for a brief review, see Kaaret & Ford 
1997). The existence of the QPOs can be exploited 
to measure additional properties of the x-ray emission 
mechanisms and geometry. Here, we use kHz QPOs 
with relatively high coherence, u/Au > 100, to mea- 
sure microsecond lags between emission in different 
x-ray energy bands. 

In §2, we describe the observations and technique 
used to calculate the time lags. In §3, we present the 
results of the lag measurements. We conclude, in §4, 
with brief comments on the physical implications. 

2. Observations and Analysis 

The following results use data from the Propor- 
tional Counter Array (Zhang et al. 1993) on the Rossi 
X-ray Timing Explorer (Bradt et al. 1993) with 122 [is 
time resolution for x-rays in the 2-70 keV band. The 
data presented here for 4U 1636-536 are from a 23 ks 
interval beginning UTC 1998 Feb 24 23:26:20 during 
which QPOs were detected. 

Before attempting to measure time delays, we in- 
vestigated the QPOs. We performed Fourier trans- 
forms on 1 s segments of data with no energy selec- 
tion, yielding a Nyquist cutoff frequency of 4096 Hz 
and a transform window function with a width of ap- 
proximately 1 Hz. The power spectra within 128 s 
intervals were summed. The total power spectrum 
for each 128 s interval was searched for QPO peaks 
above 100 Hz. A function, consisting of a Lorentzian 
plus a constant, was fitted for each trial frequency, de- 
fined as those points with power more than 2a above 
the average. All 128 s intervals for which QPO peaks 
were detected with greater than 2a significance were 
retained. In the case of multiple QPO detections in 
one interval, the QPO with highest significance was 
chosen. Typically, the QPOs were detected at better 
than 4a significance. 

The single QPO peak had a centroid varying in 
the range of 775-895 Hz, see figure |l]. The centroid 
frequency could be determined to an accuracy better 
than 0.4 Hz. The QPO peak widths had an average 
of 4.7 Hz. The intrinsic variation in the centroid is 
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Fig. 1. — Centroid of the QPO peak versus time for 
4U 1636-536. The time zero is UTC 1998 Feb 24 
23:26:20. 

far larger than the centroid measurement error or the 
QPO peak width. 

In using an oscillation (whether coherent or quasiperi- 
odic) to measure a time delay between different en- 
ergy bands, the cross spectrum is integrated over a 
selected frequency range. A narrow frequency range 
reduces the noise in the measurement, but the range 
must be sufficiently broad to contain the oscillation 
power. In the standard approach, a fixed frequency 
range is used for the full duration of an observation 
(e.g. Vaughan et al. 1997). Due to the large vari- 
ations in the QPO centroid in our observation, we 
chose instead to integrate the cross spectrum over a 
range varying with the QPO centroid. We selected an 
integration range of approximately twice the average 
QPO peak width. Compared to the fixed frequency 
range of 120 Hz required to contain the wandering 
QPO peak, the signal to noise ratio for measurement 
of a time delay in our observation of 4U 1636-536 is 
improved by roughly a factor of 3. 

We calculated the cross spectrum between a pair 
of energy bands for each 1 s of data and found the 
average and variance for each 128 s interval. A com- 
plex cross correlation vector was calculated for each 
128 s interval by integrating over a 10 Hz frequency 
range centered on the QPO peak centroid found for 
that interval. We summed the complex cross cor- 
relation vectors found for each interval to derive a 
complex cross correlation vector for the whole ob- 
servation. The phase lag given by the argument of 
this total complex cross correlation vector was then 
converted to a time delay using the average of the 
QPO centroid frequencies found in the 128 s inter- 
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Fig. 2. — Time delay versus energy for 4U 1636-536 
(squares) and 4U 1608-52 (diamonds). The error bar 
on energy indicates the band used to calculated the 
delay. The delays are relative to the 3.8-6.4 keV 
band for 4U 1636-536 and to the 3.9-6.3 keV band 
for 4U 1608-52. Following Cui et al. (1998), the neg- 
ative values emphasize the fact that the hard X-rays 
actually lead the soft X-rays. 

vals. The error was calculated from the variance in 
the real and imaginary components of the total cross 
correlation vector. Intrinsic as well as statistical fluc- 
tuations contribute to the error estimate. We also es- 
timated the time delay for the whole observation by 
averaging time delays calculated for each 128 s inter- 
val. The two methods gave consistent results, which is 
expected since the fractional range of the QPO cen- 
troid variation is relatively small, about 14%. We 
cannot distinguish whether the lags remains constant 
in phase or in time as the centroid frequency varies. 
Below, we quote the time delay found from the total 
cross correlation vector. 

To determine the correct sign for the lags, we tested 
our analysis code using simulated data and observa- 
tions of SAX J1808. 4-3658 for which the correct sign 
can be determined by direct examination of the pulse- 
folded light curve in the various energy bands (Cui, 
Morgan, & Titarchuk 1998). 

3. Results 

Time delays, for 4U 1636-536, were calculated from 
138 segments of 128 s each in which QPOs were de- 
tected. The soft emission in the 3.8-6.4 keV band 
lags the hard emission above 9.3 keV by 25.0 ± 3.3 [is. 
The coherence, calculated following the prescription 
of Vaughan & Nowak (1997), is consistent with unity 



and is above 0.75. The detection of the lag is signif- 
icant at the 7cr level. The lag for the energy bands 
6.4-9.3-17.0-29.4 keV versus the 3.8-6.4 keV band 
are shown in figure |[ We indicate the lags as nega- 
tive to emphasize that the hard emission precedes the 
soft emission. The magnitude of the delay increases 
with energy. 

We also divided the data into the four contiguous 
segments of QPO detections apparent in figure [l[ We 
find no strong evidence for variability of the time de- 
lay between any two segments. The most significant 
difference is only 1.6a. 

To check our analysis and also to have a second 
source for comparison, we analyzed data for 4U 1608- 
52 for a 13 ks interval beginning UTC 1996 Mar 
3 19:18:20 extracted from the public archive. The 
4U 1608-52 data show a single QPO varying over the 
range 820-890 Hz. We find the same sign and similar 
magnitudes of delay in 4U 1608-52 as in 4U 1636-536. 
We detected QPOs in 66 segments of 128 s each and 
used an integration width of 12 Hz to better match 
the QPO width in 4U 1608-52. The soft emission in 
the 3.9-6.3 keV band lags the hard emission above 
9.5 keV by 21.4 ± 2.7 [is. The energy dependence of 
the lag is shown in figure [|. Delays are shown for the 
energy bands 6.3-9.5-17.4-29.0 keV relative to the 
3.9-6.3 keV band. The energy bands were selected to 
best match those used in the 4U 1636-536 analysis, 
allowing for gain changes in the PCA. Vaughan et al. 
(1997) found delays of similar magnitude, but oppo- 
site sign for part of this same observation. The sign 
was later corrected (Vaughan et al. 1998) and is now 
consistent with our results. 

4. Discussion 

The facts that the hard emission precedes the soft 
emission and that the modulation of the QPO is 
higher in the hard band (Zhang et al. 1996) suggest 
that the process which generates the QPO acts ear- 
lier and more strongly at higher energies. The soft lag 
could be produced either by reprocessing of the hard 
emission into a softer band or by spectral softening of 
the intrinsic emission during each oscillation cycle. 

The lag is opposite in sign to that generally ex- 
pected in Comptonization models, in which the spec- 
trum is usually generated by upscattering soft pho- 
tons (Sunyaev & Titarchuk 1980), but could be due 
to Compton scattering of hard photons injected into 
a scattering medium which is cooler than the injected 
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photon spectrum. This mechanism has been sug- 
gested for the x-ray millisecond pulsar, SAX J1808.4- 
3658 (Cui et al. 1998). Using the delay between the 
3.8-6.4 keV and 9.3-17.0 keV bands leads to an es- 
timate (Cui et al. 1998) of ~ O.lrkm for the size of 
the Comptonization region with optical depth r. For 
r ~ 12, as typically found for 4U 1636-536 (Christian 
& Swank 1996), the Comptonization region would be 
very compact, at most a few kilometers. 

Additional observations are required to measure 
more accurately the variation of the lags with energy 
in order to determine if the soft lags are generated 
by Comptonization or by intrinsic softening of the 
emission on each QPO cycle. Also more realistic the- 
oretical modeling, including generation of the input 
hard photons, is required to adequately interpret the 
observations. 

P. K. and S. P. acknowledge support from NASA 
grants NAG5-7405 and NAG5-7334. 
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